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SERVOMECHANISMS,  EXPONENTIAL  SMOOTHING,  AND 
A  MULTIECHELON  INVENTORY  PROBLEM 


1.  Servomechanisms 

1.1  Review  of  Previous  Results 

In  a  recent  report  [1]  the  problem  of  designing  an  inventory- 
system  from  a  servomechanism  point  of  view  was  examined.  The  inventory 
model  adopted  was  a  periodic  review  system  with  fixed  order  points, 
taken  as  the  nonnegative  integers  for  the  sake  of  convenience.  The  t 
demand  period  is  measured  by  real  time  between  t-1  and  t  for 
t  =  1,2,3, •••  •  Adopting  a  fixed  time  lag  T  for  delivery,  an  order 
is  to  be  placed  at  the  end  of  each  demand  period  in  anticipation  of 
future  demands. 

It  is  natural  under  the  above  circumstances,  to  suppose  that  the 
til 

demand  during  the  t  demand  period  is  described  by  a  discrete  parameter 

stochastic  process  X.,  t  =  0,1,2,...  ,  where  X  =  0.  Furthermore, 

for  t  >  1,  it  is  assumed  that  X^  =  m(t)  +  e^,  which  describes  the 

process  in  terms  of  its  mean  value  function  m(t)  and  an  error  term. 

The  errors,  {e^},  are  considered  to  be  uncorrelated  -with  E[e^]  =  0 

and  V[e^]  s  a  .  If  we  define  the  inventory  level  over  a  safety  level 

th 

at  the  end  of  the  t  period  to  be  1^,  then,  allowing  back  orders 

against  future  deliveries,  a  simple  recursion  relation  exists  relating 
% 

inventories  to  demands  and  orders,  namely, 

(1.1)  =  Xt-1  +  ®t-l-T  '  Xt  t-  1,2,3,... 
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In  this  relationship, 


I  =0  and  0  =0  for  m  <  0.  It  is  seen  that, 
o  m  — 


admitting  no  control  over  demand,  the  inventory  level  is  automatically 
controlled  once  an  ordering  rule  0^  is  specified.  It  is  this  automatic 
control  feature,  employing  the  feedback  of  previous  inventory  levels, 
that  naturally  suggests  an  analysis  from  a  servomechanism  point  of  viev. 

Adopting  as  a  minimum  requirement  that  an  order  rule  should  be 
chosen  in  such  a  vay  as  to  control  the  inventory  about  a  safety  level, 
order  rules  vere  further  restricted  to  those  vhich  could  be  represented 
as  linear  combinations  of  past  demands  and  inventory  levels.  Thus,  for 
some  sequences  (A^A^Ag, . . . }  and  {BQ,B1,B2, . . . },  ve  require 

a'2>  *.  ■  i  Vw  *  I,  Vw  “W-- 

Next,  recognizing  that  if  is  random  then  so  is  1^,  the 
concept  of  minimizing  inventory  level  vas  replaced  by  that  of  minimizing 
the  expected  value  of  the  inventory  level  denoted  =  E[I^],  t  =  0,1,2,...  . 
Because  of  the  amenability  of  the  form  of  equations  (l.l)  and  (1.2)  to 
z -transforms,  equivalent  expressions  vere  then  derived  to  yield, 
respectively, 

(1.3)  l(z)  =  zl(z)  +  zT+10(z)  -  X(z)  and 

(1.4)  0(z)  =  A(z)X(z)  +  B(z)l(z)  . 

In  the  above  notation,  if  (V0,V^,Vg,  •  •  • }  is  a  sequence,  its  z-transf orm 
00  ^ 

V  V.z^  vill  alvays  be  denoted  V(z).  Solving  (1.3)  and  (1.4) 
k=0  K 
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r 


simultaneously,  it  is  possible  to  express  the  basic  relationship  between 
inventory  and  demand  as, 


(1.5) 


l(z)  a  S(z)X(z),  where  S(z) 


zT+1A(z)-1 

1-z-zT+1B(z) 


S(z)  is  often  called  the  transfer  function  and,  if  preferred,  a  similar  re 
lationship  exists  between  the  expected  inventory  level  and  the  mean 
demand,  namely. 


•  (1.6) 


E(z)  =  S(z)m(z)  . 


S(z)  is  the  same  transfer  function  as  defined  in  (1.5)* 

The  problem  is  thus  reduced  to  determining  an  order  rule  by  means 
of  specifying  the  sequences  {A^A^Ag, . . . }  and  {BQ,B1,Bg, . . .  }  or, 
equivelently,  their  transforms  A(z)  and  B(z).  Now  for  m(z)  =  1, 
i.e. , 


if  t  =  0 
if  t  >  0 


i 


would  be  the  response  in  expected  inventory  to  a  one-time  impulse 
at  t  =  0.  In  this  case,  E(z)  =  S(z),  and  the  requirement  that  the 
response  eventually  dampen,  i.e.,  lim  E .  =  0  or  lim  S.  =  0  is 

tw  .  *C 

->  00  t  ->  CO 

equivalent  in  the  transform  domain  to  the  requirement  lim  (l-z)s(z)  =  0. 

z  -» l” 

Since  this  in  turn  will  be  true  if  the  poles  of  S(z)  are  outside  the 

_  4*r 

unit  circle  in  the  z -plane,  we  follow  Vassian  [2]  and  let 
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•Vf; 

i 


l-z 


In  that  case,  S(z)  has  no  finite  poles  and,  with  this  determination 
of  B(z),  all  that  remains  is  to  choose  A(z).  Mathematically  it  turns 
out  to  be  easier  to  select  a  certain  function  of  A(z)  namely 


A*(z)  =  zT+1A(z)  X(z)  . 


After  some  algebraic  simplification  and  inversion  from  the  transform 
domain,  it  is  possible  to  express  (A£)  as  the  sequence. 


(!.7)  At+l+r  It+1-+T 


t+l+T 


-  £  V 

j=t+i  ° 


t  =  0,1,2, .. .  and 


A*  =  0  for  j  =  0,1,2,..., T 


Taking  expectations. 


(1-8) 


EtAt+i+r1 


t+l-KT 

Vi+t  +  J=E+1  m(J>  • 


In  view  of  the  basic  role  played  by  A^  in  (l. 7 ),  it  is  quite  suitable 

that  A*+l4T  be  called  the  forecast  of  demand  from  period  t+1  through 

period  t+l-HT.  Moreover,  the  requirement  that  lim  E.  =  0  forces 

t  Z 

t+l+T 

as  t  -» ».  This  condition  is  referred  to  as  asymptotically  accurate 
forecasting. 


A*+i+t  to  be  chosen  in  such  a  way  that  E(A*+l4^)  -  £  m(j)  -*  0 

J=tH 


| 


Finally,  once  an  asymptotically  accurate  forecaster  is  selected, 
ve  can  easily  express  the  order  rule  in  terms  of  such  a  forecast  by 
means  of  the  relation, 

<1-»)  ®t  ■  AW+I  -  I,  9t-J  -'t  *>-T  ' 

«K 

In  terms  of  servomechanism  diagrams,  the  inventory  system  can  be 
visualized  quite  easily  as  in  Figures  1.1  and  1.2  below,  first  as  a 
simple  input-output  system  and,  secondly,  as  a  more  complicated  system 
with  a  feedback  loop  clearly  displayed. 


Figure  1.1:  Input -Output  System 


Figure  1.2:  The  Inventory  System 
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1.2  Additional  Results  for  the  Case  of  Constant  Mean  Demand 


Whenever  the  demand  process  is  assumed  to  have  a  constant  mean, 
it  is  possible  to  obtain  several  explicit  results  vhich  in  turn  provide 
further  insight  into  the  inventory  system  just  described.  To  this  end, 
ve  assume  henceforth  that  m(t)  =  a  where  a  >  0  so  that,  in  trans¬ 
form  notation,  m(z)  =  .  Given  the  demand  history  X^  ...,X^ 

up  to  time  t  >  T,  ve  first  use  exponential  smoothing  to  estimate  mean 
demand  thereby  obtaining 


(1.10) 


vhere 


p  =  1  -  a,  and  0  <  a  <  1 


The  constant  a  is  called  the  smoothing  constant.  It  is  then  natural 
to  define  the  forecast  of  demand  over  the  following  T+l  time  periods  by 


At+HT  =  (T+1)at 


thereby  obtaining 

EtA?+l+T]  =  (T+1)EKJ  =  8(T+l)(l-P'°)  . 


In  this  case, 

t+l+T  . 

E^At+l+T^  ’  ^i  =  -a(^+l)P  -»0  as  t 

as  required  for  asymptotic  accuracy. 

Having  selected  an  appropriate  forecast  one  can  then  find  an 
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-***+*£&  r^XT' 


explicit  formula  for  A(z).  Indeed,  since 

A 

a(z)  = 


QX(z) 
1-P  z 


it  follows  that 


A.(.)  •  (M)»^») 


But  also. 


A*(z) 


z^U-z1*1) 

ITS 


A(z)X(z) 


so  that,  upon  substitution. 


(1.11) 


A(z) 


g(T+l)(l-z) 

(i-pz)(l-zT+1) 


As  a  casual  observation,  with  the  above  explicit  formula  for  A(z),  we 
can  express  the  transfer  function  explicitly  as 


1  PZ  k=0 


Since 


E(z)  =  S(z)m(z)  =  , 


it  follows  that 

lim  _  (l-z)E(z)  =  a  lim  _  s(z)  =  a[(T+l)-(T+l)]  =  0 
z  — ►  l”  z  -» l" 


as  required. 
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With  the  present  form  for  the  transfer  function  it  is  also  possible 
to  express  the  inventory  level  strictly  in  terms  of  past  demands.  Since 


/  v  T+l  T+l 

I(z)  =  S(z)X(z)  =  X(z)  -  X(z) 


it  follows  from  inversion  that 


(1.12) 


I.  =  a(T+l) 


t-T-1 

Jb 


Pkx 


t-T-l-k 


-  I  * 


k=0 


t-k 


for  t  >  T 


Also, 

4  A 

(1.13)  et  =  a(T+l)[Xt-at_1J  +  Xt  for  t  >  1  . 

The  above  formulas  will  be  useful  in  the  ensuing  discussion. 

1.3  Variance  of  Inventory  Level 

While  the  requirement  of  asymptotic  accuracy  has  a  great  deal 
of  intuitive  appeal  and  represents  a  minimum  requirement  in  adopting 
a  forecasting  technique,  it  is,  after  all,  a  mean  value  property.  Even 
though  the  expected  inventory  level  is  effectively  zero  after  a  certain 
time,  it  does  not  of  course  follow  that  the  actual  inventory  level  is 
even  near  zero  at  that  time.  This  fact  was  dramatically  demonstrated 
recently  by  DeWinter  [3]*  Simulating- the  above  mode],  random  demands  were 
generated  and  the  corresponding  random  inventory  levels  were  observed 
for  many  time  periods. 

In  one  example  treated  by  DeWinter,  the  expected  inventory  level 
was  effectively  zero  after  83  time  periods.  This  was  a  case  where 
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a  =  16  and  T  =  4.  Exponential  smoothing,  with  a  constant  of  a  =  0.2 
vas  used.  The  system  was  allowed  to  run  for  1,000  time  periods  using  a 
random  input  corresponding  to  a  Poisson  demand  with  mean  16.  Even  after 
900  periods,  the  actual  inventory  level  was,  at  times,  as  low  as  -33  units 
and  as  high  as  +27  units.  Other  examples  treated  by  DeWinter  showed  a 
similar  behavior  in  the  inventory  fluctuations. 

The  above  results  are  in  a  sense  not  too  surprising  since  the 
variance  of  the  inventory  level  had  never  been  considered  explicitly. 

With  the  formula  (1.12)  derived  above,  it  is  now  possible  to  analyze 
the  inventory  variance  expl.citly.  Recall  that  in  the  present  model 
we  have  assumed  X^  =  m(t)  +  where  the  random  variables  {e^}  are 
uncorrelated  with  constant  variance.  In  that  case  (1.12)  may  be  used  to 
compute  the  variance,  V^,  of  1^  as  follows.  From  (1.12)  we  may 
write 


where,  for  t  >  T 


r  a(T+l)pt-T-1"J 

L  -i 


for  0<j<t-T-l 
for  t  -  T  <  J  <  t  . 


Since  X  =  0, 
o 


and  consequently, 

V.  =  <72  T  7?  =  ^(T+l)2  Y  p2(t-T-1"«3)  +  (t+1)ct2  . 

*  J=1  J  A 
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After  some  simplification, 


(1.1*0  Vt  =  a2(T+l)[l  +  (l-p2(t-T“1))]  . 

Now,  for  positive  K,  the  function 

-  m 

defined  for  0  <  x  <  1  has  a  derivative  equal  to 


2  2K 

Since  x  <  1  and  x  <  1,  it  follows  that  f*(x)  <0  so  that  f 

is  monotone  decreasing  in  x.  This  means  that  is  monotone  decreasing 

2 

in  P  or  monotone  increasing  in  a.  Moreover,  for  a  -  0,  =  cr  (T+l) 

2 

while  for  a  =  1,  =  <r  (T+l)(T+2).  Consequently  we  have  found  finite 

bounds  for  the  inventory  variance,  valid  for  all  t  >  T 
* 

(1.15)  ct2(T+1)  <  Vt  <  <j2(T+1)(T+2)  . 

Allowing  t  to  increase  beyond  bound  in  ( 1. 14)  the  limiting  variance, 
V,  of  the  inventory  level  is  found  to  be 

(1.16)  V  =  cr2(T+l)[l  +  • 

Of  course  V  is  monotone  decreasing  in  p  as  is  V^.  As  a  matter  of 
curiosity,  taking  the  constants  used  by  DeWinter  in  the  above  example, 
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i.e.,  T  =  4,  a  =  .2  and  a  =  16,  we  find  a  limiting  variance  of 
V  =  125*  Now  three  standard  deviations  would  be  roughly  33  units. 
Consequently  for  a  value  of  t  sufficiently  large,  an  inventory  level 
of  -33  (c.f.  p.  8)  would  be  compatible  with  the  asymptotic  variance. 

The  results  of  the  simulations  are  thus  quite  in  keeping  with  the  theory. 
Moreover,  since  =  E[I^]  =  -a^+l)^  <  0  we  would  expect  the  actual 
inventory  levels  to  tend  toward  the  negative  side  of  the  zero  mean.  Such 
appears  to  be  the  case,  not  only  In  this  example,  but  in  others  examined. 


2.  Exponential  Smoothing 


2.1  Biasedness 

For  the  simple  problem  described  in  Section  1.2,  the  exponentially 

A 

smoothed  estimate  of  a  is  given  by  equation  (1.10).  The  estimate 
can  be  represented  recursively  by, 

(2.1)  2t  =  £Wt+pat_1  (t  =  1,2, ... )  , 

A  A 

which  exhibits  a^  as  a  linear  combination  of  a^_1  and  X^  with  the 
present  observation  receiving  "weight"  a. 

A 

Note  from  (2.1)  that  computation  of  a^  does  not  require  knowledge 

A 

of  the  entire  demand  history  but  only  a^_1  and  X^.  As  a  result,  if 
computations  are  performed  on  an  electronic  computer,  memory  is  conserved. 
This  feature  along  with  the  intuitive  appeal  of  weighting  the  present 
more  than  the  past  and  the  computational  simplicity  accounts  for  a  large 
part  of  the  popularity  of  the  procedure  as  an  estimating  technique. 

It  is  interesting  to  note  however  that 

E[at]  =  (l-St)a  for  t  =  1,2,... 

/N 

That  is  to  say,  the  exponentially  smoothed  estimator  a^.  underestimates 
a  for  all  finite  t,  although  it  is  unbiased  in  the  limit.  The  amount 
of  bias  is  greatest  for  small  t  and  in  particular 

A 

E[a1]  =  Cta  . 

When  a  is  chosen  between  .1  and  .3,  as  is  often  suggested,  the 
initial  error  can  be  substantial. 


Based  on  the  first  observation,  our  forecast  of  demand  during  the 

o 

next  T+l  periods  is 

■  (T+1K  < 

which,  on  the  average,  will  underestimate  future  demand.  The  resulting 
order  would  then  be  inadequate  to  satisfy  demands  occuring  during 
the  lead  time  and  an  unfavorable  inventory  position  will  occur. 

Figure  2.1  shows  1^  for  a  simulated  facility.  Demand  at  the 
facility  in  period  t  was  generated  according  to  a  Poission  distribution 
with  m(t)  s  16.  Lead  time  T  was  assumed  to  be  4  periods.  Using 
procedures  described  in  Zehna  [1], 

At+1«  -  <T+1>“  t  e\-k  (t  •  0,1,2, ... ) 

k=0 

and 

0t  =  At+HT  "  ^  0t-j  It  ^  • 

Figure  2.1  shows  that,  after  starting  with  Io  =  0,  1^  falls 

sharply  until  period  6  when  the  first  order  is  received.  Unfortunately, 
the  biased  character  of  the  forcast  resulted  in  this  order  being 
inadequate. 

Subsequent  orders  are  similarly  inadequate  and  recovery  from  an 
unfavorable  position  (-80  units)  is  seen  to  be  slower  than  would  be 


desirable. 


2.2  Finite  Exponential  Smoothing 

In  order  to  remove  the  bias  inherent  in  exponential  smoothing, 
a  modification  is  proposed. 

Let  X^,Xg, .  ..,X^  denote  a  sequence  of  observations  from  a  process 
having  mean  value  function  m(t)  =  a.  Define  the  finite  exponentially 
smoothed  estimate  of  a  by 


(2.2) 


vhere 


*  “t  I. 


X  a  o,  a.  =  -r  ,  0  <  a  <  1  ,  and  p  =  1  -  a  . 

0  *  i-p*  -  ~ 


It  follows  that  has  the  desirable  property  that  for  all  t 


E[at]  =  a  . 


That  is,  Sj.  is  unbiased  and,  a  fortiori,  asymptotically  unbiased. 
Analagous  to  expression  (2.1)  is  the  recursion  relation 


(2-3) 


at  =  at  [Xt  +  0/at.1)St.1)  (t  =  1,2,...) 


which  follows  from  the  expression  (2.2).  Since  1  +  =  l/a. , 

^  A 

equation  (2.3)  shows  that  a^,  like  a^,,  is  a  weighted  average  of  the 

A 

present  demand  observation  and  the  past  estimate  of  a.  As  with  a^, 
a.  is  simple  to  compute  and  does  not  require  that  the  demand  history 

It  will  be  convenient  to  abbreviate  "exponential  smoothing"  by 
e.s.  and  "finite  exponential  smoothing"  by  f.e.s. 
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remain  in  computer  memory  if  computations  are  computerized. 
From  the  definitions  of  e.s.  and  f.e.s.  it  follows  immediately  that 


(2.4) 


lim  [a.  ]  =  lim  E[a.  ]  =  a 
t  — » «  t  -*  » 


2.3  An  Order  Rule  for  Finite  Exponential  Smoothing 

One  of  the  significant  consequences  of  the  automatic  control 

feature  of  the  model  presently  being  examined  is  the  fact  that  once  a 

forecaster  is  given,  the  order  rule  is  automatically  determined  and  so 

then  is  the  inventory  level.  This  is  not  to  say,  however,  that  an 

explicit  formula  for  6 ^  or  1^  is  easy  to  obtain. 

We  have  just  seen  that  if  finite  smoothing  is  used  for  estimating 

the  constant  *ea„,  a,  then  the  est^tor  a,  is  unbiased  for  all  t. 

However  it  should  be  noted  that  in  the  form  at  =  st 

th  k“0 

is  expressed  as  the  ordinary  product  of  the  t  term  of  the  two  sequences 

{a^J  and  In  this  form,  an  explicit  expression  for  the 

k=0 

z-transform,  a(z),  is  not  easy  to  obtain.  In  fact,  the  z-transform  of 
CL 

a  = - -  iS  itself  not  a  recognizable  one  in  closed  form.  Also,  while 

1-P 

it  is  still  quite  natural  to  forecast  by  means  of  the  formula 


(2.5)  ^t+l+T  =  (T+1)  ®t  9 

there  is  no  convenient  or  even  recognizable  form  for  the  z-transform, 
A*(z).  Looking  over  the  development  of  Section  1.2,  it  is  apparent 
that,  without  such  an  explicit  formula  for  A*(z),  neither  A(z)  nor 
S(z)  can  be  expected  in  closed  form  by  the  same  technique.  In  turn, 
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this  prevents  an  explicit  formula  for  the  inventory  level  1^  from  being 
derived  by  means  of  inverting  transforms. 

It  is  however  possible  to  obtain  an  order  rule  analagous  to  (1.13) 
directly  from  the  recursive  relations. 

Letting 

(2-6)  9t  -  *?+1+T  +  £  \  t  =  , 

so  that  =  <p. ,  it  can  be  shown  by  induction  that 

(2.7)  dt  =  q>t  -  for  t  >  2  . 

Upon  substituting  expressions  (2.3),  (2.5)  and  (2.6)  into  (2.7)  we  obtain, 
for  t  >  1, 

(2.8)  9t  =  ^t(T+l)(Xt-at_1)  +  Xt  , 

the  analogue  of  expression  (1.13).  Since  a^  and  a^  are  asymptotically 
equivalent,  it  follows  that  0^  and  0^  are  similarly  related.  Hence, 
for  large  values  of  t,  the  order  quantities,  based  on  equivalent  fore¬ 
casting  techniques,  will  be  the  same  as  with  e.s.,  and  consequently, 
the  inventory  levels  will  also  be  the  same. 

2.4  Use  of  Finite  Exponential  Smoothing 

In  Figure  2.2  the  performance  of  the  simulated  facility  when 
using  e.s.  is  compared  with  that  using  f.e.s.  The  simulation  operated 
identically  except  for  the  different  smoothing  procedures.  The  random 
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sequence  of  demands  vas  the  same  for  both  problems  in  order  to  facilitate 
a  comparison  of  results. 

/N 

If  we  denote  by  1^  the  inventory  at  the  end  of  period  t  vhen 
e.s.  is  employed  as  a  forecast  proceudre,  and  by  the  corresponding 
variable  vhen  f.e.s.  is  employed,  ve  note  that 

1^  =  ^  ( t  =  1,2, . . . ,5 ) 

At  the  end  of  period  6,  hovever,  our  inventory  position  is  sub¬ 
stantially  better  using  f.e.s.  This  improvement  is  a  result  of  the 
unbiased  character  of  the  f.e.s.  forecast.  Periods  1  through  24  vividly 
show  the  rapid  recovery  obtained  with  f.e.s. 

There  are  suggestions  in  the  literature  that  vhen  a  change  in  the 
process  generating  the  data  is  suspected  one  should  disregard  data  which 
originated  prior  to  the  change,  since  it  is  no  longer  pertinent  to  fore¬ 
casting  the  future.  If  this  is  done  one  must  employ  f.e.s.  to  avoid  the 
bias  that  would  result  from  starting  anew  and  using  e.s.  as  the  forecast 
procedure. 

Figure  2.2  also  shows  that  for  periods  176  through  200, 


This  result  is  expected  for  t  sufficiently  large  in  view  of  expression 
(2.4). 


The  inventory  at  the  end  of  period 


(2-9) 


t-T-1 


t  is 

t 
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Substituting  expression  (2.7)  into  (2.9)  and  simplifying  the  summation 


we  note  that 


h  =  '  J&  ^ 


Substitution  for  9  and  simplification  yields 


(2.10) 


I.  = 


t-T-1  .  T 

(T+l)at_T_1  Y.  p  Xt-T-l-k“  E  Xt-k 
k=0  k"^) 


Equation  (2.10)  is  the  f.e.s.  analogue  of  equation  (1.12). 

Following  the  procedure  described  earlier  for  1^,  we  can  obtain 
^  the  variance  of  1^.  Thus, 


’*  ■  &  'A 


where 


r  (T+1)«t-T-i  p 


t-T-l-j 


l<J<t-T-l 
t  -  T  <  J  <  t 


Hence 


~2 

vt 


[(T+l)  +  (T+l)2o^_T_1  t^1  p£(t"T'1'J)]a2 


which  upon  simplification  yields 


2(t-T-l) 

(2.11)  -  cr2(T+l)[l  +  (T+l)(_|±_5)(i±_ - )]  . 

(1-P  / 


For  K  >  1,  the  function 
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i 


I 


. 


-  (gic)(^) 


defined  for  0  <  x  <  1  has  derivative 


f'(x)  = 


-2  +  2KxK~1  -  2KxK+1  +  2x2K 
(l-xK)2(l+x)2 


To  show  that  f*(x)  <  0  it  is  sufficient  to  show 


l-x2K  >  KxK"1(1-x2) 


or,  equivalently. 


Noting  that 


1  +  x  +  •••  +  x2K-1  >  KxK~1(  1+x )  . 

x  is  a  convex  function  of  K,  it  follows  that,  for 


1  <  r  <  K, 


x 


r-1 


xr-l  <  xK+(K-r) 
+  x2K'r  >  x^1 


+  x 


xr-l+(K-r) 


or 


Summing  over  r. 


Y  xr  1  >  KxK’^(l+x)  as  required  . 
r=l 


This  follows  from  the  fact  that  if  f  is  convex,  y1  <  y2  and 
h  >  0,  f(y2)  -  f(y1)  <  f(y2+h)  -  f^+h). 
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AVrn  . 


W 


^2 

Thus,  as  would  be  expected,  is  monotonically  decreasing  in  p. 

For  p  =  0,  =  cr2(T+l)(T+2)  while,  for  P  =  1,  L* Hospitals  rule  is 

employed  to  yield  =  a2(T+l)[l  +  (T+l)(^~^-)].  Hence  for  all  t  >  T, 
~2 

has  finite  bounds 

(2.12)  <t2(T+1)(1  +  (T+l)(^j)l  <  ?2  <  ^2(T+l)(T+2)  • 

Allowing  t  to  approach  infinity  in  expression  (2.1l)  yields  the 
limiting  variance 

(2.13)  V  =  <72(T+1)[1  +  . 

/>2  A 

Denoting  by  the  variance  of  1^  and  recalling  that 


f 


it  follows  that  for  t  >  T 

(2.14)  V2  >  . 


Letting  t  increase  without  bound, 


lim 
t  -*  00 


t 


02 

lim  V2 
->  00 


Expression  (2.14)  would  seem  to  be  in  conflict  with  the  simulation 
results,  which  indicated  that  the  fluctuations  about  zero  were  smaller 
for  f.e.s.  than  e.s.  The  findings  are  however  compatible  since  the 
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deviation  about  zero  for 
is,  if  Et  =  E[It], 


includes  the  effect  of  a  bias  term. 


That 


A  p  A  A  A  p 

E[Ij7  ]  =  E[It-Et+Et]c 


=  V 
vt 


/\2 
+  Et 


while 


3.  Multiechelon  Supply  System 

3.1  A  Servomechanism  Model  of  the  System 

The  system  just  described  may  be  generalized  to  a  multiechelon 
system  by  casting  it  into  the  framework  of  a  multiple  input -output 
servomechanism.  As  an  example,  consider  a  system  composed  of  tvo  tenders 
supplied  from  a  common  depot.  Here  there  are  two  inputs,  namely,  the 
demands  xj.1^  and  X^2^  from  the  two  tenders  and  there  are  three  outputs 
of  interest;  these  are  the  respective  inventory  levels  and  l£2^ 

at  the  tenders  along  with  the  inventory  level  1^^  at  the  depot.  Of 
course,  the  demand  at  the  depot  level  is  the  sum  of  the  two  order  rules 
and  initiated  at  the  two  tenders. 

Considering  each  of  the  tenders  to  be  operating  according  to  the 
model  outlined  in  the  previous  section  each  of  the  tender  inventory  levels 
may  be  expressed  by  means  of  z -transforms  in  terms  of  separate  transfer 
functions.  Thus, 

(3.1)  I^(z)  =  S^(z)X^(z)  and 

I^(z)  =  S^2'\z)X^2^(z)  . 

Now  thinking  of  the  depot  as  constituting  still  another  version  of  the 
same  model  with  an  input  of  0^  +  there  is  a  transfer  function 

S^(z)  such  that 

(3.2)  i^\z)  =  s^(z)(e^(z)  +  e^2hz))  . 


Because  of  the  relationship 


.  i 


e^\z)  =  A^UJX^U)  +  B^(z)l^(z)  i  =  1,2 

it  is  possible  to  express  I^\z)  further  in  terms  of  the  demands  at 
the  tenders.  In  fact,  since 

I^(z)  =  S^(z)X^(z)  , 


it  follows  that  for  i  =  1,2 

e^\z)  *  A^^zjx^^z)  +  B^^zjs^^zjx^^z) 
=  [A^(z)  +  B^^zJS^^z)]  X^(z)  . 


This  result  allows  us  to  express  all  of  the  input-output  relationships 
simultaneously  in  matrix  form.  Following  the  notation  of  Howard  [4], 
we  obtain. 


(3-3) 


where 


— 

I(l)(s) 

hn(z) 

0 

I(2)(z) 

= 

0 

^22^  ^ ^ 

I(3)(z) 

h5l(z) 

h52(z) 

— m  — 

_ 

(z) 


hn(z)  =  S(l)(z)  , 

h2g( z )  *  S^S\z)  , 

h51(z)  =  S^\z)[A^(z)  +  B^UJs^U)]  and 
h32(z)  =  S(5)(z)[A(2)(z)  +  B(2)(z)S(2)(z)]  . 
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For  visual  aid,  a  flow  chart  of  this  multiechelon  system  is  given 
in  Figure  3.1.  We  have  assumed  that  T^  is  a  delivery  lag  time 
i  =  1,2,3,  and  in  general  T^,  Tg  and  T^  are  distinct.  Observe  that 
the  order  rule  at  the  depot  level  0^\ z),  necessary  to  maintain  its 
own  inventory  level,  is  a  linear  combination  of  its  input  (the  sum  of 
the  tender  orders)  and  its  invenotry  level.  In  fact, 

Q^\z)  =  A^^(z)[e^1\z)  +  Q^2\z)]  +  ^\z)l^\z) 
as  in  the  single  echelong  model. 

3*2  Reduction  to  Single  Echelon  Model 

The  final  remark  made  in  the  previous  section  suggests  the 
possibility  of  analyzing  the  inventory  at  the  depot  level  as  simply 
another  single  echelon  model.  This  can  in  fact  be  accomplished  but  not 
without  a  little  difficulty.  Assume  that  the  demands  at  the  tenders 
have  constant  means  so  that  there  exist  constants  a1  and  a 2  such 
that  a.^  =  m^t)  =  EtX^J  for  i  =  1,2,  whence  m^^z)  =  and 
( 2 )  a2 

mv  '(z)  =  are  the  respective  transforms  of  the  mean  value  functions. 
Assume  further  that  exponential  smoothing  is  used  to  forecast  demands 
at  each  of  the  facilities.  Without  loss  of  generality  we  suppose  that 
each  facility  is  required  to  use  the  same  smoothing  constant  a. 

Now  the  results  of  Section  1  apply  to  each  of  the  tender  levels, 
and  we  obtain  as  in  that  section. 
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where 


Rk  Jl) 
@  Xt-k 


is  the  smoothed  estimate  of  the  constant  mean  a^  From  this  it  follows 

% 

immediately  that 

(3-5)  Efe^]  =  ai  +  a1Q{(Ti+l)pt"1  i  =  1,2. 

Consequently,  E[0^  +  9 =  (a^ag)  +  a3t"1[a1(T1+l)  +  a2(T2+l)], 

a  quantity  which  depends  upon  t.  This  is  unfortunate  since  +  0^ 

is  the  input  at  the  depot  level.  Thus,  even  though  the  simple  model 
of  constant  mean  demand  is  imposed  at  the  tender  levels,  it  does  not 
follow  that  the  same  model  applies  at  the  depot  level.  Because  of  the 
nature  of  the  input  at  the  depot,  namely,  the  expected  value  being 
an  exponential  function  of  time,  the  analysis  in  terms  of  z-transforms 
is  not  particularly  suitable. 

On  the  other  hand,  if  finite  exponential  smoothing  were  used  at 

the  tender  level  then  the  expected  orders  would  be  constant  and  the 

th 

simple  model  would  then  apply.  Accordingly  we  require  that  the  i 
tender  place  orders  according  to  the  rule, 

+  40 


where  aj,*| 
a  given  by 


is  the  finite  exponentially  smoothed  estimate  of  mean  demand 

si-i -  °t-i  ki  pM-L-  As  tefore’  at =  •  in 


that  case,  we  have 
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f  kJ  . . si 


El^1^  =  ^(Tj+lJIa^a^  +  8i  =  8i  , 


i  =  1,2 


so  that 

E(®11)  +  ^  -  al  +  a2  ’ 
a  constant  independent  of  time. 

Letting  x[^  =  denote  the  input  to  the  depot,  the 

t  t  t 

latter  may  now  operate  as  a  single  echelon  model  with  constant  mean 

demand  The  depot  order  rule  6^^  will  then  be  determined 

(3) 

once  the  sequence  A'  or,  alternatively,  the  depot  demand  forecast 

(3)* 

At+l+T  is  determined.  In  this  regard,  the  depot  now  has  a  choice 

between  finite  and  ordinary  exponential  smoothing  as  a  forecasting 

technique  and  it  should  be  observed  that  in  either  case  the  resulting 

(3) 

order  rule  '  is  a  function  of  second  order  exponential  smoothing 
applied  to  the  tender  demands  and  xj2^.  This  is  because  the  depot 

input  X^  is  already  a  weighted  sum  of  tender  demands. 

t 

The  investigation  of  the  variance  of  the  depot  inventory  level 
is  not  so  simple.  For  one  thing  the  single  echelon  model  does  not  quite 
fit  in  this  case  for,  even  though  the  expected  demand  is  constant,  the 
variance  is  not.  In  fact,  if  f.e.s.  is  used  at  the  depot  level  then, 
as  previously  derived, 

=  at(Ti+1^xt1^“®t-i^  +  Xt^  for  1  =  lf2 


or. 


8[1)  -  [“t<v1)+llxt1)  -  at(Ti+1)at-i  £  pkxt-i- 

k=0 
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J'K'i— 


in  which  case. 


V[8^]  =  [a  (T^lJ+llV  +  °^at_i(Ti+1)2  £  ^ 

k=0 


2t 

=  a2[Q^(Ti+l)2  +  2at(Tj+l)  +  1  +  Q;2a2_1(T1+l)2  g?.]  for  i  =  1,2. 

1-P 

Now  it  is  quite  reasonable  to  suppose  that  the  order  rules  ad) 

and 

*(2) 

0^  from  the  two  tenders  are  independent.  In  that  case, 

V[X^]  =  V[s[1'1]  +  V[9^2)]  =  o^to^Tj+l  f  +  2at(T1+l)  +  Q^(Tg+l)  +  2at(Tg+l) 
+  2  +  fa&p  £<£  ((T  +1)2  +  (T  +1)2)] 

d-e2)  tW  1  2 

which  is  not  a  constant  although  lim  V[xP^3  ia  constant.  What  is 


t  — »  » 

worse,  since  both  2^^  and  $£2^  are  linear  combinations  of  their 

respective  demand  histories  up  to  time  t  it  is  not  evident  that  the 

(3 ) 

members  of  the  sequence  (X£  '}  are  uncorrelated.  Hence,  the  computation 
of  the  variance  of  the  inventory  level 


I<3)  =  a(T  +1)  *  J'1  ekx[-T  1-k 
k=0  5 


-  I 


k=0 


,0) 

t-k 


(for  t  >  T^)  is  not  at  all  straightforward.  Thus,  simulation  techniques 
will  have  to  be  adopted  to  estimate  the  variance  of  the  depot  inventory 
level. 
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3  •  3  Generalization 

The  simple  3  facility  system  just  described  can  be  generalized 
to  one  composed  of  n  facilities.  Material  is  shipped  between 
facilities  according  to  a  supply  procedure  specified  by  a  set  of  integers 
fpfg, ...,fn  with  f^  =  0  and  1  <  f j  <  J  for  j  >  1.  The  integer 
f j  is  interpreted  as  the  unique  facility  directly  supplying  facility 
J.  Let  F.  =  {k|f,  =  J}  (j  =  l,2,...,n)  be  the  set  of  facilities 
directly  supplied  by  facility  J. 

It  is  convenient  to  assume  that  the  basic  time  unit  is  the  period 
and  that  all  units  of  time  can  be  expressed  as  an  integral  number  of 
periods . 

0 

Each  facility,  j,  employs  a  periodic  review  policy  whereby  at  the 
end  of  Rj  time  periods  it  evaluates  its  stock  position  and  if  appro¬ 
priate  places  an  order  with  facility  fj.  The  order  is  delivered  to 

facility  J,  T.  periods  later. 

~( l) 

Let  X£‘  denote  the  user  demand  at  facility  J  in  period  t,  and 


let 


v  _  rvC1)  $(n)\ 

At  ”  '  At  9  9  *  *  * 9  ' 


be  the  vector  of  those  demands.  If  a  particular  facility  has  no  user 
demand  but  only  functions  to  supply  other  facilities  then  ^  =  0. 

In  general  we  assume  that  the  user  demand  process  can  be  expressed 


J  =  1, 2, . . .  f  n 


where 


mj(t)  is  the  deterministic  mean  value  function 
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and 


e^(t)  is  a  random  variable  with  E[£j(t)]  =  0. 


Let  J  =  { j | Fj  -  0 } .  Then  for  j  £  J  let 

x(j).,_  j(j) 

t  ■  t  ^ 

and  for  j  0  J. 

(3.6)  X(J)  =  3j(J)+  £  £  e“  j=l,2,...,n 

z  t  meF,  u+K  £(t-R,+l.t) 
j  m  J 

where 

0™  represents  an  order  from  facility  m  made  at  the  end  of  period 
u, 
and 

represents  the  lag  between  when  an  order  is  made  at  facility 

m  and  the  stock  is  removed  at  facility  f  . 

m 


Thus  if  Fj  ^  0  the  demand  imposed  on  facility  J  is  composed  of  the 
user  demand  and  the  demands  which  have  arrived  within  the  review  period 
from  the  facilities  which  it  supplies.  In  many  supply  systems  expression 
(3*6)  will  simplify. 

The  previous  example  is  a  special  case  of  this  formulation  in  which 


xt1)  s  *t2-'  42M3) 

f2  =  f3  =  ! 

F1  =  (2,3),  F2  =  0,  F5  =  0  , 
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J  =  {2,3} 


f 


3*4  Depot  Forecasting  Information 

The  simulation  described  earlier  was  extended  to  the  above 

~(2) 

three  facility,  tvo  echelon  supply  system.  It  was  assumed  that 
and  *?>  were  Poisson  distributed  random  variables  with  parameters 
\2  =  \^  =  1,  and  »  Tg  *  Tj  =  4.  In  accordance  with  Section  3.2, 

(3-7)  *t+l  *  ^  +  ^  (t  >  1)  .. 


Simulation  was  used  to  examine  the  relationship  between  system 
performance  and  the  nature  of  the  information  transmitted  from  facilities 
2  and  3  to  facility  1. 


Case  I: 

In  this  case  we  assume  that  the  only  Information  being  received 
by  facility  1  is  the  order  information.  In  particular,  we  assume  that 
a  forecast  is  obtained  from  expression  (2.3),  where 


*U) 

t+1 


is  given  by  expression  (3*7). 


Case  II: 

As  an  alternative  to  this  forecast  procedure  at  facility  1,  we 
examined  the  situation  in  which 

a*1)  =  3t£2)  +  , 

ill  -  + ill ♦  (*  >  o) 
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In  Case  II,  facility  1  makes  its  forecast  from  the  demand  information 
originating  at  facilities  2  and  3.  We  assume  that  this  information  is 
supplied  to  facility  1  without  delay. 

Demand  information  supplied  without  delay  has  two  advantages.  It 
provides  a  preview  of  the  order,  to  be  imposed  on  facility  j 

during  the  next  period  and  thereby  effectively  reduces  T.  It  provides 
demand  information  whereas  0^  describes  demand  information  "confounded" 
with  inventory  adjustments. 

Because  of  the  above  advantages  we  might  expect  better  performance 
in  a  system  which  supplies  demand  information  to  its  supplier  as  it 
originates. 

Two  hundred  periods  of  system  operation  were  simulated.  The 
performance  of  facility  1  as  measured  by  the  end  of  period  inventory 
is  shown  in  Figures  3*2  and  3*3  for  periods  1  through  50  and  150  through 
200  respectively. 

Figures  3*2  and  3*3  clearly  indicate  reduced  variation  in  Case  II. 
Further  investigation,  both  analytical  and  with  simulation  techniques 
is  desirable  to  examine  the  nature  and  source  of  the  variation. 
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“  Period  Inventory  for  Facility  1  in  a 
Facility  Multiechelon  Supply  System 


End  ol'  Period  Inventory  for  Facility  1  in 
3  Facility  Multiechelon  Supply  System 
Periods  151-200 


I 


& 


4.  Conclusions 

4.1  Servomechanism  Techniques 

Servomechanism  techniques  continue  to  be  helpful,  not  only  in 
finding  explicit  formulas  for  such  quantities  as  order  rules,  but  also 
for  obtaining  a  much  better  understanding  of  the  inventory  system  and 
the  manner  in  which  it  operates.  Analyzing  the  model  treated  in  this 
paper  from  such  a  point  of  viev  has  served  well  to  discover  certain 
weaknesses  as  well  as  strengths  in  the  system.  Certainly  it  has  greatly 
emphasized  the  important  role  played  by  forecasting  techniques.  This 
type  of  analysis  has  also  led  to  recognition  of  the  trade-off  in  exponential 
smoothing  between  response  to  changes  in  demand  patterns  and  the  desire 
to  minimize  variation  in  inventory  levels.  Also,  some  progress  has  been 
made  in  the  direction  of  structuring  a  multi-echelon  system  using  vector¬ 
valued  functions.  The  reduction  of  such  a  complex  system  to  a  more  trac¬ 
table  one  is  another  reason  why  this  approach  should  be  exploited  in  the 
study  of  inventory  systems. 

4.2  Finite  Exponential  Smoothing 

The  analytical  results  and  the  confirming  simulation  indicate 
the  virtues  of  finite  exponential  smoothing  as  a  forecast  technique, 
a^  was  seen  to  be  unbiased  for  all  t.  This  meant  a  better  inventory 
position  during  the  first  several  periods  following  any  initialization  of 
the  demand  history. 

a^  can  be  computed  recursively  with  two  relations,  thus  requiring 

/\ 

only  little  more  effort  or  computer  storage  than  a^. 
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5.  Further  Research 

ts 

2  2  ^ 

We  know  that  by  making  £  large  a  (X^),  c  (i^)  and  consequently 
cyclic  variation  can  be  reduced.  The  penalty  in  so  doing  is  however 
the  loss  of  sensitivity.  This  occurs  since  the  exponentially  smoothed 

A  A 

forecase  is  a  linear  combination  of  X^.  and  with  1-0  as 

the  weight  given  to  the  current  observation.  Therefore,  if  0  is  large 

A 

little  change  occurs  in  a^  with  X^. 

Consider  a  simple  problem  in  which  the  mean  function  m(t)  of  the 
demand  process  {X( t ) )  has  the  following  form 


m(t)  =  a  1  <  t  <  t 

o 

m(t)=a  +  c  t>t 

—  o 


In  a  situation  such  as  this  it  would  be  desirable  to  employ  f.e.s.  with 
0  close  to  1  throwing  away  all  data  occuring  prior  to  the  change. 

Such  a  procedure  would  require  a  mechanism  for  detecting  changes 
in  the  demand  pattern.  Such  procedures  have  been  suggested  in  the 
literature  [5],  [6]-  The  application,  evaluation,  and  modification  of 
such  procedures  to  the  forecasting  problem  represents  an  important  area 


for  further  research. 
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